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ABSTRACT 

We report on the N II AA5668-5712 emission and absorption lines in the 
spectrum of rj Carinae. Spectral lines of the stellar wind regions can be classified 
into four physically distinct categories: 1) low-excitation emission such as H I 
and Fe II, 2) higher excitation He I features, 3) the N II lines discussed in this 
paper, and 4) He II emission. These categories have different combinations of 
radial velocity behavior, excitation processes, and dependences on the secondary 
star. The N II lines are the only known features that originate in "normal" 
undisturbed zones of the primary wind but depend primarily on the location of 
the hot secondary star. N II probably excludes some proposed models, such as 
those where He I lines originate in the secondary star's wind or in an accretion 
disk. 

Subject headings: stars: emission-line, Be - stars: individual (eta Carinae) - stars: 
variables: general - stars: winds, outflows 

1. Introduction 

In i] Car's broad-line stellar wind spectrum, the high-excitation helium features have dif- 
ferent profiles and fluctuate differently from the "normal" lines of H I, Fe II, etc. Most authors 
now assume that the H e I em ission depends on photoionization by a hot companion star, see 



6 of [Humphreys et al.l (120081 ) and references therein. Observed He II emission is most likely 
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excit e d by soft X-rays from un stable shocks (Martin et al.ll2006l . but cf. ISteiner fe Damineli 
2004 . Kashi fe Soker 2007 . and Soker fe Behai 12006 ). Features of these types are important 
because they show direct influences by the secondary star and the wind-wind collision region. 
The only clear examples have been helium lines, whose source geometry is both complex and 
controversial. In this paper we report similar characteristics in a set of N II features, which 
sample lower-ionization gas than the helium lines. 

Broad H I and Fe II emission lines represent p Car's stellar win d spectrum that would be 
present even if the companion star did not exist (iHillier et al.ll2001l). Their main components 



alway s remain close to system velocity (roughly —8 km s 1 , 
20041 ). 



Davidson et al. 1997 



Smith 



The He I emission lines, by contrast, shift progressively blueward through most of the 
5.5-year spectroscopic cycle. When a spectroscopic event occurs they show a more abrupt 
negative shift, followed by a rapid positive reversal to renew the cycle (see mid panel of 
Figure HJ and figures in iNielsen et al.l 120071 ). The overall range of variation is about to 
—250 km s _1 . Evidently these changes represent flows of highly ionized material, modulated 
by the secondary star. As the secondary moves in its orbit, it progressively illuminates 



regions with differing velo city fields, and 



dependent: compare , e.g., iNielsen et al 



Kashi fc Soker 2007, 2008 
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them. The details are model- 
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2001 



a. 



He I absorption shows a similar pattern with variations between —300 km s~ x and 
—550 km s^ 1 . This behavio r is mirrored by the puzzling He II A4687 e mission observed 
during spectroscopic events (ISteiner fe Daminelil 120041 ; iMartin et al.l I2006T ) . Hydrogen and 
Fe II lines have components that behave qualitatively like He I, but cannot be measured 
separately because they overlap the steady "normal" emission profiles^ 

Until now no other species were found in r\ Car with velocity variations like the helium 
lines, and without conspicuous steady components. Here we report that lines of the N II 
AA5668-5712 multiplet do match this prescription. They are always present but were not 
discussed previously because they were weak. They are noteworthy because they depend on 
a form of excitation by the hot secondary star, but in different regions than the He I lines. 
They seem reasonable in some models, but are difficult to explain in others. For example, 
models where the helium emission originates in a small region close to the secondary star 
have difficulties in this respect; see §4 below. 



1 This is particularly clear in STIS data showing US (Fig. 5 in lNielsen et al.l 120071 ). Velocities quoted 
above a re unpublished meas ure ments of Gemini GMOS ob servations during the 2009 spectroscopic event, 
see also lNielsen et al.l (|2007l ) and lSteiner fe Daminelil (|2004l ) for similar values during the 2003.5 event. 
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Historically, hot windy objects such as P Cygn i show strong N II emission and absorption 



Beals 1935; Struve 



1935 



Swings k, Struyd Il940f ) . Since nitrogen is abundant in rj 



(e.g., 

Car's CNO-processed wind, iHillier et al.l ( 120011 ) expressed surprise that N II is not mor e 
conspicuous there. Some of the f eatures discussed below were identified by iGaviolal ( 119531 ). 
bu t were too weak to be noted by iThackeravl (119531). More recently they have been detected 
by lZethsonl ( 1200ll ). iDamineli et al.l ( 119981 ). and lNielsen et al.l (120091 ) without discussion. The 
lines discussed here are physically different from [N II] A5756 which has much lower energy 
levels. 



2. Data and Analysis 



Ground-based slit spectroscopy of rj Car was obtained with the Gemini Multi-Object 
Spectrograph (GMOS) on the Gemini South telescope from 2007 June to 2010 January. Here 
we only use 20 sets of observations with the B1200 line grating and a 0.5" slit that cover 
the spectrum at A5700 A on the star and on FOS4@ FOS4 is a position in the SE lobe of 
the Homunculus (about 4" southeast of the star) which reflects the star's polar spectrum 



while the spectrum in direct view is representative of lower la titudes (IDavidson et al.lll995 



Zethson et al.lll999t IDavidson et al.ll2001bl ; ISmith et al.l 120031 ). The spectral resolution was 



about 1.2 A or 65 km s -1 . The seeing varied from 0.5" to 1.5", so each GMOS spectrum 
represents a region typically ~ 1" across. We prepared 2-D spectrograms with the standard 
GMOS data reduction pipeline in the Gemini IRAF package and extracted 1-D spectra. The 
pipeline wavelength calibration was improved using the interstellar absorption line at A5782 
A; in this paper we use heliocentric vacuum velocities. Further details, unimportant for this 
paper, will be listed in a later publication (Mehner et al. 2011, in preparation). 

Eta Car was observed with the Hubble Space Telescope Space Telescope Imaging Spec- 
trograph (HST STIS) from 1998 to 2004 and then again starting in 2009. Here we only use 
observations with the G750M grating, the 52" x 0.1" slit, and central wavelength A5734 A. 
The spectral resolution was about 1.1 A or 60 km s _1 . The 1998-2004 observations include 
a variety of slit positions and orientations. 2009 June and December data cover a region up 
to 1" from the central source with parallel slit positions offsets of 0.1". The long-exposure 
2009 data were obtained in HST programs 11506 and 12013, whose Pi's were K. S. Noll and 
M. Corcoran. The data were reduced using routines developed by the HST Treasury Project 
that include several improvements over the normal STScI pipeline and standard CALSTIS 



2 For a list of these observations and other details see the HST Treasury Program website 
p:/ /etacar.umn.edu/| in particular Technical memo number 14. 
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reductionsEI We extracted spectra 0.1" wide. 



3. Results 

In the GMOS data we found broad emission and absorption lines of the N II 2s 2 2p3s 
3 P° - 2s 2 2p3p 3 D multiplet at AA5668-5712 A, exhibiting radial velocity variations during 
the 2009 spectroscopic event similar to the helium lines; see Figure [1] and Table [TJ Figure 
[2] shows the energy levels involved. Note that the lower level, 2p3s 3 P°, is more than 18 eV 
above the N II ground level but is connected to it by a permitted transition. We discuss the 
excitation mechanism below. 

Figure [1] compares spectra of r] Car obtained from 2007 June to 2010 January showing 
N II A5668, He I A4714, and He II A4687. Weak N II AA5668-5712 emission and absorption 
is always present in the spectrum of rj Car. The maximum strength of the absorption 
component has a radial velocity of about —300 km s _1 with respect to the emission peak. 
During the 2009 spectroscopic event, emission and absorption features shifted about 250 km 
s _1 towards the blue, simultaneous with the He I emission and absorption lines and the He II 
A4687 emission. 

HST STIS data obtained in 2009 June and December clearly show the N II AA5668-5712 
lines, see Figure [3j In retrospect the faint N II spectral features can also be detected in HST 
STIS data from 1998-2004, but due to lower S/N they failed to attract notice before. Like 
the He I absorpti on, N II absorption w as much weaker before 2003; less than 30% compared 
to 2009 (compare iMehner et al.ll2010bl ). The absorption strengths of both species increased 
gradually over the last 10 years except during the 2003.5 and 2009.0 spectroscopic events. 
We determined the "stellar continuum" by a loess (locally weighted scatterplot smoothing) 
fit to the spectrum and find that th e N II absorp tion is stronger than the emission. (The 
same is true for P Cygni; lBealslll935l . IStruvdll935l ). Gemini GMOS data appear to indicate 
the opposite, but their data quality is too low to determine a reliable continuum, and ground- 
based spectra are contaminated by narrow emission lines from ejecta outside the wind, e.g., 
[Fe II] A5675 blends with the N II AA5678,5681 absorption. Although we first noticed these 
features in GMOS data, the high spatial resolution of HST data is essential for examining 
their character. 



Most other permitted N II lines are too weak or blended with emission lines from the 
wind or nearby ejecta to be observed. However, we see similar behavior in the N II 2p3s 3 P° 



For information see http://etacar.umn.edu and iDavidsonl (2006) 
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- 2p3p 3 P series at AA4603,4608 A (Fig. [2]). Transitions of the N II 2p3s 3 P° - 2p3p 3 S series 
are too weak and blended with other lines. 

Radiative excitation of level 2s 2 2p3s is very important as discussed later. The 2s2p 3 
levels are excited in the same way but lead to no observable results. Transitions from 
2s2p 3 3 S° to 2s 2 2p3p (AA6435-7265 A) have small oscillator strengths and are therefore not 
observed. STIS data indicates the weak presence of the singlet N II A3996 line but this line 
is blended with others in Gemini data and therefore the identification is not certain. Lines 
whose upper levels are above the 2p3p state like N II A5007, even though they are strong in 
the laboratory, are not observed. 

The N II lines and their velocity shifts can also be seen in the reflected polar spectrum 
at location FOS4, the location decribed in §2. The velocity shifts from the FOS4 separate 
direction may cast doubt on orbit models, but this question is too complex to discuss here. 



4. Line Formation and Significance 



The same N II lines are seen in early type stars. In objects like P Cygn i and shells of O- 
stars, 2p3s 1 P° - 2p3p 1 S, 1 P, 1 D and 2p3s 3 P° - 2p3p 3 S, 3 P, 3 D are very strong (jSwings fc Struve 
19401 ) . Large differences between WN stars indicate that these lines are sensitive to atmo- 
spheric conditions and/or the variability of the wind. They may also be very sensitive to the 
FUV flux since they are likely produced via continuum fluorescence (IHerald et al.ll200ll ). 



In the case of rj Car the N II absorption and emission almost certainly occur in the 
primary star's dense wind, for reasons noted later below. But three facts indicate that UV 
photons from the secondary star populate the critical 2p3s level. (1) The velocity behavior 
described in §3 strongly suggests some relation to the secondary, analogous to the He I 
features. (2) The 2p3s levels are abou t 18.5 eV ab o ve th e ground state, a high value for 
the primary star's wind. According to iHillier et al.l (j200lf ). the opaque- wind photospheric 
temperature is below 15, 000 K and emission-line regions in the primary wind are mostly 
below 10, 000 K; much cooler than the O stars and WR objects mentioned above. The 
ho t secondary star, on the other hand, produces a large flux of 18.5 eV photons (Fig. 10 



m iMehner et al.l l2010al ) . (3) This hypothesis appears quantitatively successful as outlined 
below. Any of these clues by itself might be debatable, but together they seem compelling 
to us. 

Let us summarize an order-of-magnitude assessment of the absorption line strengths that 
one would expect in a simple model. For simplicity we include only the N + ground level '1' 
(2p 2 3 P) and one excited level '2' (2p3s 3 P°); our initial goal is to estimate the equilibrium 
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population ratio n^jny. A two-level system is justified because no permitted transitions 
connect level 2 to intermediate levels. Consider a uniformly-expanding local volume in the 
primary wind. Denote the incident continuum photon flux at hv^. ~ 18.5 eV by not 



corre cted for line absorpt ion in the gas. As explained in §8.5 in lLamers fc Cassinelhlll999l and 



8 in iMartin et al.l 12006k the Sobolev approximation allows us to write expressions for the 



radiative 1 — > 2 excitation rate and the re-emission photon escape probability, as functions 
of the local expansion rate, a line profile function, and atomic parameters. The effective de- 
excitation rate is proportional to the escape probability, and the most complicated factors 
appear similarly in both the excitation rate and the escape probability. Consequently the 
equilibrium population ratio is simple: 



U2 = #2 _ A^ 2 $ 

ni 



X 



9i 



71 



where g\ and #2 are the statistical weights. This expression remains valid when we take fine 
structure into account. The same ratio would be found in an LTE case where <& v /c is the 
Planckian photon density at /iz/ 12 rj 18.5 eV. Above we mentioned triplet levels of N II, but 
singlet levels would also be excited in the same way from 2p 2 X D. 

With conventional parameters for the two stars and their winds, <& u at 18.5 eV is dom- 
inated by the hot seco ndary star. Reasonab le values for it are T e g rs 40, 000 K, R « 13 R® 
and L rs 4 x 10 5 (IMehner et al.ll2010al ). According to a WM-basic atmosphere model 



(IPauldrach et al.ll200ll ). such an object emits roughly 5 x 10 33 photons Hz 1 s 1 at hv = 18.5 



eV, about 30% less than a 40, 000 K blackbody. Figure H] shows th e probable geometry . 



With the type of orbit model that most authors currently favor (e.g., lOkazaki et al.ll2008l ). 
the separation between stars was about 13 AU when rj Car was observed in 2009 June. The 
secondary star was then located roughly 1-5 AU closer to us than the primary, and roughly 
10 AU from our line of sight to the primary - depending on the poorly known orbit orienta- 
tion, of course. With these parameters, eqn. 1 gives n%jn\ ~3x 10~ 8 at relevant locations 
between us and the primary star, i.e., in gas located about 10 AU from the secondary star. 
The equivalent excitation temperature is near 12, 500 K. 

The observed absorption lines have some geometrical complications. Since the opaque- 
wind photosphere is diffuse with a substantial size, a relevant "line of sight" can be any 
sample ray in a cylinder with diameter ~ 7 AU (Fig. HI). Strong 2p3s-2p3p absorption 
occurs in regions that optimize a combination of attributes: (1) the inner wind is favored 
because of its high densities; (2) most of the nitrogen must be singly ionized; (3) the N + 
must be fairly close to the secondary star in order to maximize <&„; and (4) the light path 
(possibly indirect, due to Thomson scattering) must pass close enough to the primary star 
to sample continuum radiation from its semi-opaque inner wind. Thus, in Fig. HJ we expect 
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the strongest absorption to occur near and above the symbol '1' that labels Region 1. The 
picture obviously changes as the secondary star moves along its orbit. 

Most nitrogen in the primary wind is s ingly ionized, due to both the primary and the 
secondary radiation field (IHillier et al.ll200ll ). Helium in zones 3 and 4 of Fig. H] provides the 
shielding which prevents the secondary radiation field from ionizing N + into N ++ , but which 
allows the 18.5 eV radiation to penetrate. Suppose the primary mass-loss rate is 3 x 10 M Q 
yr _1 and nearly all the CNO is nitrogen ( Davidson et al. 1986 ; Dufour et al. 1999b Fl Then 
the column densities outside r = 4 AU (for example) are 



iV(N" 



2 x 10 20 cm" 2 . 



6 x 10 12 cm" 2 



N(N 2p3s) 

This column density would produce a combined equivalent width of about 1 A for the 
A5681.14 and A5677.60 A absorption lines. Since this exceeds the observed value of 0.5 A, 
the proposed mechanism appears to be sufficient, even if only a limited part of the wind is 
involved. Let us emphasize that the parameter values in this sample calculation were chosen 
ab initio without knowing what result they would lead to; there was no readjustment to get 
a desired outcome. Note that the margin of error does not allow much smaller values of $> u 
and gas densities. 

The emission features reported in §3 are comparable in strength to the absorption - 
probably weaker and certainly not much stronger - so they resemble the pure scattering 
case of P Cyg lines. In other words, the observed 2p3s - 2p3p emission line strengths 
automatically resemble the absorption strengths. More detailed calculations will require 
elaborate geometrical models of the ionization and excitation zones. 

As we noted earlier, these N II features should arise chiefly in regions of the primary 
wind that are fairly close to the secondary star, and, therefore, close to the He + zones and 
wind- wind shocks. If shock instabilities do not distort it too much, the He + /He° ionization 
front is expected to have a pseudo-hyperboloidal shape that may be either concave or convex 
towards the secondary star (Fig.HJ). Adjoining primary- wind zones are spatially large enough 
to account for the N II lines. Generically, this type of model can explain the velocity 



variations, but authors disagree about deta ils (e.g.. iNielsen et al.ll2007t iDamineli et al.ll2008 
Humphreys et allhood : iMartin et aDl2006h . 



4 The most often quoted mass loss rate for ij Car is around 10 3 M P) yr 1 , but there are good reasons 
to believe that it has decreased substantially in the past decade ( Mehner et al. 2010b ; Corcoran et al. 20ld 



Kashi fc Sokerll2009h . 
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Far more important than merely being consistent with some models of the rj Car sys- 
tem, the N II features may disprove others. For example, it has been suggested that the 
helium lines originate in the acceleration zone of the secondary star's wind, or perhaps in an 



accretion disk around the secondary, rather tha n in the primary wind ( ISoker fc Beharl 12006 



Kashi fc Sokerl 120071 ; ISteiner fc Daminelill2004j ). But the N II features, with practically the 



same velocity behavior, almost certainly cannot represent such zones; so models of that type 
would be forced to explain the He I and N II velocities differently. 

The arguments are simples t for absorption lines. First we note tha t with credible pa- 



rameters for the two stars (e.g., iHillier et al.ll200ll ; iMehner et al.ll2010al ). the primary wind 



accounts for at least 95% of the visual-wavelength continuum, most likely 98-99%. There- 
fore the absorption features in Fig. definitely represent material located between us and the 
inner parts of the primary wind. (In other words, blocking the entire secondary star would 
not produce absorption as deep as that shown in Fig. [3j even if rii ne ^> 1.) Can this absorbing 
material be part of the secondary wind or an accretion disk of the secondary star? Presum- 
ably not; because in every proposed orbit model the relevant lines of sight either miss the 
secondary star and its wind entirely at most times, or pass through only the outskirts of the 
secondary wind. Gas densities in the secondary wind are two orders of magnitude less than 
we used in the calculations outlined above, the nitrogen is mainly N III, and the velocity 
dispersion is too large for the observed absorption features. This type of model would badly 
fail the quantitative feasibility test that the "conventional" model passed. Regarding the 
N II emission features, they are fairly weak relative to the total observed continuum, but 
would be extraordinarily strong compared to the visual- wavelength continuum of just the 
secondary star. 

Developing these arguments further would be beyond the scope of this paper; our main 
point is that the N II lines convey information that is not available in other spectral features. 
This is mainly because they arise in "normal" parts of the primary wind but they are 
sensitive to the current location of the secondary star. Helium lines also depend on UV from 
the secondary, but their emission zones are extremely difficult to model since they depend 
on the unstable primary shock, post-shock cooling, local "clumping," etc. In principle the 
N II zones are expected to be much steadier, and the excitation mechanism is relatively 
insensitive to inhomogeneities. For these reasons the N II features may be good indicators of 
average density in the wind combined with the secondary star's orbital motion. In particular 
one might expect their velocities to be easier to model than those of helium lines. 

If the "conventional" model is more or less valid, then the N II absorption features 
should weaken around apastron when <& v is smallest. Unfortunately the 2000-2002 HST 
STIS data (the only previous apastron passage when STIS was available) are not adequate 
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to test this "prediction" easily. Perhaps new STIS observations will be made in 2011-2012. 

In §3 we mentioned that He I and N II absorption lines have grown in strength during 
the past 10 years. At first sight this appe ars counter-intuitive, since the wind density appears 



to have decreased (jMehner et al.ll2010bl ). However, the situation has various dependencies 
that tend to oppose each other so intuition may be a poor guide. Decreasing mass loss 
rate implies a smaller, denser opaque-wind photosphere; eqn. 1 becomes invalid if continuum 
absorption destroys the trapped 18.5 eV resonance photons before they escape; etc. This 
problem merits further attention. As the wind decreases it is conceivable that the primary 
star may eventually become hot enough to excite these lines in addition to the secondary. 
Whatever the solution, we expect the N II lines to provide different parameter contraints 
than the spectral features that are excited in more normal ways. 

In summary, spectral features in rj Car's wind can be assigned to these categories: 
1. The "normal" line s of hydrogen, Fe II, etc. In principle these can be used to analyze 



the primary wind ( iHillier et al.l 120011 ). 



2. He I emission and absorption. Most authors agree that these are related in some way to 
the secondary star, but the details are controversial (see references cited in §1 above). 

3. N II emission and absorption. These lines arise in normal-ionization parts of the 
primary wind, but they depend mainly on the proximity of the secondary star. In 
this sense they sample a new region of parameter space. In most recent models, the 
N II zones are much simpler than the ionized helium. We emphasize that the He I 
lines depend on photoionization, hv > 24.6 eV, whereas the N II lines depend on 
photoexcitation at hv ~ 18.5 eV. The ratio f v (> 25 eV)//j,(18.5 eV) depends strongly 



on th e secondary star's temperature (WM-basic atmosphere models, iPauldrach et al. 



200lh . 



4. He II A4687, an unsolved puzzle not discussed here. It depends on soft X-rays from 
the colliding-wind region, not UV from the secondary star. 

A good basic model for the primary wind and the binary system must account for at least 
the first three of these, throughout most of the 5.5-year cycle. (He II A4687 may depend on 
shock instabilities too complex for a "basic" model.) 



After an earlier version of this paper appeared on astro-ph, iKashi fe Sokerl ( 120111 ) have 
proposed an alternative interpretation of the N II lines. A full response would make this 
paper too long, but in our view they do not explain how gas in a very small locale near the 
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secondary star can produce substantial absorption lines in the spectrum of the much larger 
primary- wind photosphere (see discussion above). 
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Table 1. Observed N II AA5668-5712 transitions a 



Wavelength 






Transition 






E; 






(A) 














(cm -1 ) 


(cm -1 ) 




N II 5668.20 


2s 2 


2p 3s 


3 Pi 


-2s 2 


2p 3p 


3 D 2 


148940.17 


166582.45 


3.45c+07 


N II 5677.60 


2s 2 


2p 3s 


3 Po 


-2s 2 


2p 3p 


3 D X 


148908.59 


166521.69 


2.80c+07 


N II 5681.14 


2s 2 


2p 3s 


3 P 2 


-2s 2 


2p 3p 


3 D 3 


149076.52 


166678.64 


4.96c+07 


N II 5687.79 


2s 2 


2p 3s 


3 P1 


-2s 2 


2p 3p 


3 D! 


148940.17 


166521.69 


1.78c+07 


N II 5712.35 


2s 2 


2p3s 


3 P 2 


-2s 2 


2p 3p 


3 D 2 


149076.52 


166582.45 


1.17c+07 



http://physics.nist.gov/PhysRefData/ASD 
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Wavelength (A) 
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Fig. 1.— N II A5668, He I A4714, and He II A4687 in spectra of r] Car obtained with Gemini 
GMOS from 2007 June to 2010 January. Days before and after the 2009 spectroscopic event 
at MJD 54860 are displayed next to each tracing. Arrows indicate the zero velocities of the 
emission lines. Note that most of the 5.5 year cycle is represented by the top two and bottom 
two samples and that [Fe II] X5675 blends with the N II \\5678,5681 absorption. 
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Fig. 2. — The lowest triplet levels of N II. Odd-parity levels such as 2s 2 2p3s 3 P° are marked 
with small 'o's. In this paper we concentrate on the 2s 2 2p3s - 2s 2 2p3p transitions because 
they have observationally convenient wavelengths, see §3. 
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Fig. 3. — N II AA5668-5712 in HST STIS data from 2009 June. Zero velocity emission and 
absorption at —300 km s _1 are indicated with brackets. The dashed curve represents the 
stellar continuum. HST's high-spatial resolution shows the absorption much better than 
ground-based data. The strong N II AA5678,5681 absorption is blended with [Fe II] A5675 
in the Gemini spectra. 
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Fig. 4. — Schematic arrangement of a conventional model for the 2009 June observations. 
The line of sight to Earth lies in the plane of the Figure, and the secondary star's orbit is 
shown projected onto that plane. Regions 1 and 2 are undisturbed parts of the primary 
and secondary winds. The shocked wind-wind collision region is labeled 3, while region 4 
is a He II zone in the primary wind, photoionized by the secondary star. N II is abundant 
in region 1 but zones 2-4 are more highly ionized. Region 5 is within the opaque-wind 
photosphere at visual wavelengths. The area shown is roughly 25 AU across, but all details 
are simplified and idealized. 



